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ABSTRACT 


The  Physical  Optics  and  Exact  Inverse  Scattering  Solutions  of  this 
author  are  summarized.  For  the  Physical  Optics  Inverse  ocav.tering  Method, 
shown  are  computer-reconstructed  images  of  a  sphere  and  cylinder  from 
computed  synthetic  scattering  data,  as  well  as  a  sphere  from  experimentaly 
measured  data.  For  the  Exact  Inverse  Scattering  Method,  shown  are  computer- 
reconstructed  source  distributions  (currents)  of  a  half-wave  <\ipole  antenna, 
a  point -source,  and  two  point-sources  separated  by  one-half-wavelength,  from 
computed  synthetic  scattering  data. 
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SECTION  I 


INTRODUCTION 


This  report  is  a  continuation  of  the  final  report  [1]  to  the  preceding 
contract,  as  well  as  to  an  earlier  quarterly  report  [2]  to  this  contract. 

The  theoretical  results  obtained  to-date  where  reported  on  in  detail 
in  the  above  mentioned  reports,  and  will  thus  be  only  briefly  summarized 
in  this  report.  The  purpose  of  this  report  is  to  summarize  the  numerico- 
experimental  results  obtained  to-date. 

For  the  Physical  Optics  Inverse  Scattering  Method,  shown  are  computer 
reconstructed  images  of  a  sphere  and  cylinder  from  computed  synthetic 
scattering  data,  as  well  as  a  sphere  from  experimentally  measured  data. 

For  the  Exact  Inverse  Scattering  Method,  shown  are  computer  reconstructed 
source  distributions  (currents)  of  a  half-wave  dipole  antenna,  a  point  source, 
and  two  point  sources  separated  by  one-half  wavelength,  from  computed  synthetic 
scattering  data. 
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SECTION  II 


PHYSICAL  OPTICS  INVERSE  SCATTERING 


II. 1.  THEORETICAL  RESULTS 


The  basic  Physical  Optics  Inverse  Scattering  Identity  [1],  applicable 
to  complete  bandwidth  and  perspective  (aspect  angles)  information,  is 


y(X)  =  fie.  |  d3<  e 


(  d3<  e1***  F  — 

j  L  <2  . 


(i) 


where  k  s  2k,  the  range  and  phase  normalized  field-cross  section  p  is  related 
to  the  power  cross  section  or  by  pp*  s  ^5  ,  and  the  characteristic  function  y 
of  the  scatterer  is 


Y<X) 


-{  ’  ’ 
(  °  . 


x  e  v. 
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x  e  v„ 


(2) 


where  V  is  the  volume  of  the  scatterer. 
e 


If  the  information  aperture  is  incomplete  (i.e»,  only  finite  band¬ 
width  and/or  incomplete  perspective  data  are  available) ,  then  (1)  yields 
the  following  Fredholm  Convolution  Integral  Equation  of  the  First  Kind  [S] 


y(x)  *  fie  a(x)  =  fie.  j  d3<  e 

A 


/*-"■[?] 


(3) 


where  a(X)  Mk)  ,  an d  uhore  AOO  is  th®  characteristic  information  function, 
i.e. , 


A(K )  s 


1 

0 


p(*0  known 
p(x)  unknovm 


(4) 


This  integral  equation  (3)  can  be  solved  exactly  numerically  by  a  variety 
of  existing  techniques  such  as  the  matrix  methods  of  Ritz-Galc-rkin  [6],  the 
associated  Least  Square  Best  Estimate  method  [7] ,  ar /  associated  moments 
method  of  Harrington  [8] ,  the  Eigen-function  expansion  method  of  Toraldo  Di- 
Francia  [9],  leading  to  so-called  super-resolution,  and  the  k-space  method 
of  this  author  [10],  which  also  leads  to  super- resolution  (This  solution  will 
be  summarized  in  Sect.  III.l.  of  this  report,  since  similar  Fredholm  convolut'  « 
integral  equations  of  the  first  kind  arise  in  the  exact  inverse  problem) . 

Several  closed-form  solutions  of  (3)  for  apertures  of  specific  geometry 
have  been  obtained  by  Lewis  [11] . 

Since  the  unknown  function  y(X)  in  the  integral  equation  (3)  is  not  a 
completely  arbitrary  function,  but  a  characteristic  function  restricted  to  the 
form  (2),  the  following  closed- form  approximate  solutions  [12]  to  (3)  are 
obtained 


X<X) 


(5) 


R 


dx3  x3  X(X> 


(6) 


where  %{x)  is  a  three-dimensional  resolution  density  function  which  is  a 
measure  of  the  location  of  the  surface  of  the  scatterer,  and  z(x\,xy)  is  the 
geometrical  function  of  this  surface. 


II. 2.  NUMBRIOO-EXPERIMEHTAL  RESULTS 


The  solution  to  the  integral  equation  (3)  proposed  by  Lewis  [13]  was 
successfully  numerically  tested  for  a  sphere  by  this  author  in  1969.  This 
test  consisted  of  a  computer  implementation  of  a  special  case  version  of  !rhis 
solution,  applicable  only  to  scatterers  about  which  only  a  priori  knowledge 
of  cylindrical  symmetry  exists.  This  test  essentially  confirms  the  correctness 
of  the  basic  inverse  steering  identity  (1)  and  the  finite  aperture  integral 
equation  (3).  This  solution  was,  however,  not  pursued  further  because  of  its 
inherent  prect.lr.aA  limitations.  These  limitations  are  the  lack  of  generality 
nf  the  required  k-space  aperture  (i.e.,  the  required  aperture  is  impractical 
for  physically  realizable  radar  systems;  which  is  not  the  case  with  this 
author's  solution  5  and  6),  the  error  enhancement  introduced  by  the  process  of 
numerical  differentiation  of  noisy  data  (vis-a-vis  the  error  reduction  resulting 
from  the  process  of  integration  of  such  data  in  solutions  5  and  6),  and  the 
unapplicability  of  the  Fast  Fourier  Transform  (FFT)  to  this  solution  (which 
is  essential  if  large  amounts  of  data  are  to  be  processed  in  reasonable  time 
by  existing  computers,  yielding  thre-dimensionai  high-resolution  descriptions 
of  arbitrarily  shaped  scatterers  about  which  no  a  priori  knowledge  of  special 
geometry  exists). 

Solutions  (S)  and  (6)  where  computer  implemented  with  th6  aid  of  the  FFT 
for  arbitrarily  shaped  apertures,  r  alizable  with  existing  radar  systems.  lh*» 
computer  program  was  tested  with  the  exact  solution  of  Mie  for  scattering  by  a 
sphere,  with  a  variety  of  band  limited  aspect  angles  and  fractional  frequency 
bandwidths,  with  the  results  shown  in  fig.  1  and  3.  The  ripples  in  fig.  1  are 
due  to  the  very  small  three-dimr  -»ional  raster  of  163  data  points  for  z(xj,X2>; 
vis-a-vis  the  much  larger  two-dimensional  raster  of  1282  data  poins  for  x(xi,X2,0) 
in  fig.  3,  for  which  these  ripples  disapear.  A  full  bandwidth  three-dimensional 
display  of  a  reconstructed  cylinder  is  shown  in  fig.  2. 

Solution  (5)  was  also  tested  against  experimentally  measered  (in  unechoic 
chamber  by  The  General  Dynamics  Corp.,  Fort  Worth,  Texas)  scattering  data  from 
a  test  sphere;  the  results  are  shown  in  Fig.  4  and  5;  in  the  latter  figure  the 
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FIGURE  3 


TWO-DIMENSIONAL  DISPLAY  OF  RECONSTRUCTED  SPHERE  FROM  SYNTHETIC 

SCATTERING  DATA 


FIGURE  4 


TWO-DIMENSIONAL  DISPLAY  OF  RECONSTRUCTED  SPHERE  FROM  MEASURED 

SCATTERING  DATA 
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SECTION  III 


EXACT  INVERSE  SCATTERING 


III.l.  THEORETICAL  RESULTS 


The  basic  Exact  Monochromatic  Scalar  Field  Inverse  Scattering  Integral 
Equation  [14],  applicable  to  complete  perspective  information,  is 


‘ 


dv  lm  G  p  =  —  4;  dS«(G*  V«J>  -  6  VG*> 


which,  when  <}>  is  measured  in  the  far-field,  reduces  to 

Jr  _i,vx 

dv  Im  G  p  =  I  dn  e  kg  *(kg>  (8) 

v 

where  ^  is  the  range  and  phase  normalized  scattered  far-field. 

For  the  incomplete  perspective  information  aperture,  (8)  reduces  to 
the  following  Double  Fredholm  Convolution  Integral  Equation  of  the  First  Kind 


f  “ik*‘x 

a(X)  *  Im  G(X)  *  p(X)  =  I  do  e  kg  y<ks) 


The  time-domain  integral  equation  associated  with  (7)  is 


11 


Preliminary  investigations  of  the  time-domain  integral  equation  (10) 
have  been  initiated  by  Prof.  N.  Bleistein  [16], 


For  the  electromagnetic  vector  fields,  equations  (7)  through  (11)  have 
the  following  respective  analogue 


J  dV  Im  VG(xjx»)xJ(xf )  =  0(X> 

J  v 


(12) 


6  =^f(j)ds  [  VG**<nxH)  -  TC*Cn«H)  +  ius  G*(nxE)  ) 


(13) 


12 


III. 2.  NUMB  RICO-EXPERIMENTAL  RESULTS 


Figure  6  shows  the  first  iteration  synthetic  computer  reconstruction 
(computer  and  programmer  provided  by  Dr.  G.  Triccles,  General  Dynaaics  Corp., 
San  Diego,  California)  of  the  source  distribution  (current)  in  a  half -waye 
dipole  antenna,  as  per  (11).  Since  the  errors  are  reduced  by  a  factor  of 
two  per  iteration,  an  order  of  ten  iterations  should  suffice  for  most  practical 
problems.  It  should  be  noted,  however*  that  the  first  iteration  already 
yields  reconstruction  of  the  source  distribution  beyond  the  Classical  Rayleigh 
Diffraction  resolution  Halt. 

Figure  7  and  8  show  similar  results  for  the  8tn,  16^,  and  32Rd 
iterations  for  a  point  source,  and  two  point  sources  one-half  wavelengh 
apart,  respective! ly. 
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FIGURE  7 


RECONSTRUCTED  POINT  SOURCE 
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